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ABSTRACT 


An  analysis  is  developed  for  the  boundary  layer  flows,  with 
dissociated  and  ionized  air  species  in  chemical  equilibrium  or 
with  finite  rate  chemistry,  about  general  body  geometries  at 
zero  angle  of  attack.  The  partial  differential  equations  for 
the  boundary  layer,  about  axisymmetric  and  two-dimensional 
bodies,  are  transformed  to  the  von  Mises  stream  function  coor¬ 
dinate.  The  resulting  equations  are  reduced  to  an  algebraic 
system  obtained  by  adopting  an  explicit  finite-difference 
method  of  solution.  These  equations  provide  the  basis  for  a 
computer  program  to  obtain  numerical  solutions. 

Laminar  and  several  turbulent  viscosity  models  are 
included  in  the  analysis.  Turbulent  diffusivity  models  based 
upon  the  law  of  the  wall,  law  of  the  wake,  a  curve-fit  to  pipe 
flow  eddy  viscosity  and  laminar  viscosity  contributions  are 
introduced.  Constant  Prandtl,  Lewis  and  Schmidt  numbers  through 
the  thickness  of  the  boundary  layer  are  assumed.  The  thermo¬ 
dynamic  and  chemical  kinetic  data  of  the  air  species  0a  ,  0,  Ns , 
N,  NO,  N0+,  and  e~  are  considered.  Provision  is  made  in  the 
program  for  solutions  with  either  non-equilibrium,  complete 
equilibrium,  or  finite  rate  chemistry  with  the  wall  in  equi¬ 
librium. 

Boundary  layer  swallowing  of  i.nviscid  shock  layer  flew 
behind  curved  snocks  is  included.  The  streamwise  variation  in 
local  outer  edge  conditions  are  obtained  from  an  approximate 
inviscid  streamline  tracing  procedure. 
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I .  INTRODUCTION 

An  explicit  finite-difference  method  solution  has  been 
studied  for  the  equations  of  the  boundary  layer  about  hyper¬ 
sonic  reentry  vehicles  in  chemically  reactinq  air.  The 
coupling  of  the  effects  of  finite-rate  chemistry  with  various 
models  for  the  diffusional  phenomena  in  laminar  and  turbulent 


flows  was  reported  previously,  in  Ref.  1.  The  computer  program 
developed  in  the  earlier  effort  was  limited  to  the  models  of 
turbulent  viscosity,  and  only  permitted  non-equilibrium  calcu¬ 
lations.  A  somewhat  more  generalized  computational  capability 
has  been  made  available  from  the  current  investigation.  This 
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report  presents  the  results  of  this  enlargement  of  effort. 

The  capability  has  been  provided  for  solvirg  the  boundary 
layer  equations  with  the  multicomponent  air  system  in  equi¬ 
librium.  The  option  of  selecting  the  wall  in  equilibrium  and 
the  flow  in  non-equilibrium  has  been  added;  this  eliminates 
the  need  to  specify  the  wall  species  over  the  bcdy  as  input. 

An  enlarged  selection  of  turbulent  diffusivity  models  has 
also  been  included.  Two  terms  of  the  series  solutions  for  the 
analytic  determination  of  wall  parameters  and  derivatives  of 
the  functions  in  the  vicinity  of  the  wall  have  been  adopted. 

The  pertinent  boundary  layer  parameters  for  this  work 
were  defined  in  Ref.  1  and  the  methods  of  computation  were 
also  indicated  there. 
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11 •  THE  boundary  layer  equations 

The  boundary  layer  equations  for  the  steady  flow  of  a 
chemically  reacting,  multicomponent  gas  mixture,  neglecting 
radiative  energy  exchange  and  thermal  diffusion,  were  trans¬ 
formed  to  the  von  Mises  coordinates,  x-tf.  The  transformation 
from  the  x— y  plane  is  defined  bv 


e  b'i 
pur  = 

3y 


e  d# 

-  pvr  =  — 

CSX  ' 


where  the  body  radius  r(x)  is  measured  normal  to  the  center¬ 
line.  This  definition  of  the  stream  function,  %  automatically 
satisfies  the  continuity  equation.  The  body  radius  appears  for 
axisymmetric  problems  by  setting  c  =  1  and  for  two-dimensional 
problems  is  eliminated  with  c  =  C.  The  differential  operators 
become,  from  (1) , 


dx  3x 


-  Pvr 


b  e  3 

=  Pur  rr  . 
oy 


The  momentum  equation  in  the  transformed  coordinate 


becomes 


■ 
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du  _  l  dp  e  St 

dx  “  "  pu  dx  +  r  '  (2) 


0 


=  l£ 

dtf  ' 


where  the  shear  stress  is  defined  by 


T 


e  L  T 
pur  (M  '  ) 


du 

a>Irj  • 


(2a) 


(3) 


The  viscosity  coefficient,  jJ,L,T ,  may  be  either  the  laminar 
o*  turbulent  or  the  sum  of  both  in  the  computer  program 
developed  from  this  analysis. 

The  total  energy  equation  in  the  transformed  plane 
becomes 


dH 

dx 


—  (  £h 
dtf  1  c 


+  Cc) 


The  energy  f.  ,x  associated  with  the  gradients  of  total 
enthalpy  and  kinetic  energy  is  defined  by 


(4) 
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e 

pur 


dH 

d'l'  + 


e 

pur 


L,T 


(5) 
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The  Prandtl  number,  PL'T,  may  be  selected  as  any  reasonable 
value  for  the  gaseous  mixture  in  either  laminar  or  turbulent 
flows  in  the  computer  program.  The  energy  flux  associated 
with  the  diffusion  velocity  for  each  species  is  related  by 
Fick  s  law  to  the  gradient  of  the  species  mass  fractions, 
cj_»  and  species  enthalpy,  and  is  expressed  by 


f  =  pur£  1 - I_  Lj^J  (  L,T  yh  hsifh 

1  LL'Tj  \sL’Tl  '  lhi  D  9* 


where  D  represents  an  effective  average  mixture  (binary) 

diffusion  coefficient.  Laminar  or  turbulent  Lewis  numbers, 

L, T  ,  L  T 

L  ,  and  Schmidt  numbers,  S  '  ,  may  be  specified  in  the 

computer  program.  The  additional  energy  relationships  required 


are  the  static  enthalpy 


h  =  H  -  u2/2, 


and  the  static  enthalpy  in  terms  of  the  species  mass  fracti< 


and  species  enthalpies 


h  =  )  c  .  h .  . 
Z_,  i  i 


The  species  conservation  equation  in  the  transformed 


coordinates  is 


1 
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3c .  w.  .  ~ 

1  1  e  :  /  J 

3x  pu  +  r  d'l'  '  ^ 1  -  1»2#3,4...)  (7) 


where  the  net  rate  of  species  production  term,  w. ,  is  deter- 

l 

mined  by  the  phenomenological  law  of  mass  action  with  the 
experimental  reaction  rates  appropriate  to  the  production  of 
each  chemical  species.  The  term  for  the  diffusion  of  species 
is  given  by 


(8) 


A  dilute  gas  mixture  has  been  assumed  so  that  each 
species  may  be  treated  as  an  ideal  gas;  hence  the  partial 
pressure  is 


(9a) 


and  the  mixture  density  is  given  by 


(9b) 
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1 1 1  *  CHEMISTRY,  thermodynamics  and  transport  properties 

The  species  considered  in  the  (gas-phase)  air  reactions 
for  the  high  temperatures  of  interest  include  0S ,  O,  Ns ,  N, 

NO,  NO+,  and  e".  There  are  considered  active  in  the  following 
reactions  (with  catalyst  c) : 


c 

+ 

Ns 

*  2N 

+ 

C 

(10a) 

c 

+ 

Os 

-  20 

+ 

c 

(10b) 

c. 

+ 

NO 

-  N 

+ 

0  +  c 

(10c) 

NO 

+ 

0 

*  N 

+ 

0S 

(10d) 

Na 

+ 

0 

-  NO 

+ 

N 

(lOe) 

Ns 

+ 

03 

*  2  NO 

(10f) 

N 

+ 

0 

*  N0+ 

+ 

e“ 

(lOg) 

The  reaction  rates  for  the  above  reactions  are  the  same 
as  those  of  Ref.  1.  (These  laminar  rates  are  also  assumed  to 
apply  to  the  turbulent  flows.) 

Solution  of  the  equations  with  finite-rate  chemistry 
requires  the  evaluation  of  the  net  rate  of  production  of  species, 
wi'  from  the  law  of  mass  action  with  experimental  reaction  rates 
for  production  of  each  chemical  species.  The  equations  are  the 
same  as  those  of  Ref.  1.  m  order  to  analyze  equilibrium  chemical 
systems,  the  equilibrium  mass  action  laws,  rather  than  the 
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conservation  equations  (7),  must  be  solver)  for  the  species 
mass  fractions.  The  equilibrium  system  of  equations  for  the 
seven  air  species  considered  are  taken  from  Ref.  2. 

The  conservation  equations  contain  the  transport  prop¬ 
erties  of  viscosity,  thermal  conductivity,  and  diffusivity 
as  functions  of  the  flow  variables  (e.g.  temperature,  pres¬ 
sure,  and  species  mass  fractions) .  The  introduction  of  the 
nondimens ional  (constant)  transport  parameters  of  Prandtl, 

Lewis,  and  Schmidt  numbers  requires  that  only  the  viscosity 
must  be  formulated  explicitly.  The  temperature  dependency 
expressed  by  the  Sutherland  relationship  will  be  utilized  for 
the  air  in  laminar  flow.  Theoretical  descriptions  of  the 
turbulent  flow  structure  are  still  inadequate.  The  analytical 
basis  for  turbulent  boundary  layer  viscosity  models  is  dependent 
upon  semi-empirical  correlations.  Although  details  differ, 
there  is  a  general  acceptance  of  the  characterization  of 
turbulent  layers  divided  into  two  regions,  viz.,  an  inner  region 
close  to  the  wall,  and  an  outer  region  away  from  tne  wall,  with 
different  eddy  viscosity  models.  The  viscosity  model  over  the 
inner  region  of  the  flow  is  usually  taken  as  a  law  of  the  wall, 
and  over  the  outer  region  as  a  law  of  the  wake.  Several 
turbulent  viscosity  models  are  provided  within  the  computer 


program  (discussed  in  Part  II)  to  permit  combinations  to  be 


selectively  controlled.  Division  of  the  turbulent  boundary 
layer  flow  into  a  near  wall  region  with  one  viscosity  model, 
and  an  outer  region  with  another  viscosity  model,  will  be 
controlled  by  a  thickness  ratio  factor. 


KlTiMMSl'HEg!! 
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IV.  METHOD  OF  SOLUTION  -  FINITE  DIFFERENCE  EQUATIONS 

Solution  of  the  system  of  equations  discussed  in  Section 
II  by  an  explicit  finite-difference  method  was  adopted  for 
computer  programming.  The  finite-difference  formulation  was 
established  by  employing  the  explicit  difference  relations, 
for  the  characteristic  function  F,  as  follows: 


3F  _  F  (x+A  x,  ty  -F  (x,  <L') 
dx  Ax 


(Ha) 


and 


whe  re 


5F  _  F  (x,  +  A  ty-F(x,^-A 

dip  2A^ 


(lib) 


_d_  3F 


aix,^+|  M]  {f(x,  +  A^)-F(x,  “•10  } 

Ta^3 


(lie) 


a  x,  't-  y  M  {f  (x,  ty)  -F  (x,  &!jf)  } 

^—3 


a 


1X' 


=  \  [a(Xf*  ±  m  ] 


(Hd) 


The  system  of  algebraic  expressions  derived  from  these 
differential  equations  is  subjected  to  stability  conditions 
requiring  restrictions  on  the  permissible  grid  dimensions. 
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The  stability  requirements  cannot  be  precisely  formulated  for 
the  nonlinear  system  of  equations.  However,  on  the  basis  of 
linear  theory  (see,  e.g.,  Ref.  3)  and  subsequent  numerical 
testing,  the  stability  conditions  applicable  to  the  problem 
can  be  established.  The  analytical  stability  requirement  from 
linear  theory  is 

ct  A  x  .  1 

JavF  2  '  <1:?> 

where  O’  indicates  the  average  value  of  the  locally  quasi¬ 
constant  coefficients  from  the  linearized  equations.  The 
expressions  for  (X  from  the  conservation  equations  are: 


Momentum 


Energy 


Species 


O’ 

u 


2  s  .  L,  T. 
r  pu(/i  ) 


o 

H 


2  S 

r  pu 


(12a) 


(12b) 


D 


(12c) 
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V.  BOUNDARY  CONDITIONS 

The  solution  of  the  boundary  layer  equations  requires  that 
the  initial  distributions  of  velocity,  total  enthalpy  and 
species  mass  fractions  be  specified;  also,  the  boundary  con¬ 
ditions  at  the  wall  and  outer  edges  are  required.  At  the  wall 
the  boundary  conditions  may  be  stipulated  in  terms  of  the 
normal  gradients  of  the  functions,  rather  than  the  functional 
values;  in  the  case  of  the  wall  boundary  condition  of  species, 
the  element  mass  fractions  must  be  conserved.  (This  arises 
from  considerations  of  species  diffusional  velocities  at  the 
wall.)  Values  at  the  inviscid  outer  edge  of  the  layer  are 
not  known  <1  priori,  but  are  derived  from  the  boundary-layer 
swallowing  of  inviscid  flow  streamlines  (computed  within  the 
program)  as  the  calculation  progresses  in  the  streamwise 
direction.  A  streamline  tracing  procedure  for  swallowing 
behind  curved  shocks,  based  upon  a  linear  variation  of  pres¬ 
sure  with  distance  from  the  shock  to  the  body  station,  will  be 
used  to  determine  the  inviscid  flow  properties. 

For  convenience,  a  fixed  reference  coordinate  system  will 
be  adopted  to  provide  the  geometrical  descriptions  of  the  body 
contours,  shock  geometry,  relationship  of  the  body  with  respect 
to  the  shock,  and  determination  of  shock-streamline  intersection 
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points.  Various  combinations  of  spherical  sections,  ogives, 
cones,  cylinders,  flat  plates,  and  wedge  contours  will  pro¬ 
vide  the  description  of  many  body  shapes.  The  shock  s?"ape  will 
be  based  upon  combinations  of  parabolas  and  straight-line 
segments  for  curve  fits. 

The  pressure  distribution  about  the  body  is  required; 
the  source  of  such  data  may  be  experimental,  or  analytically 
determined  values  from  solution  of  the  inviscid  flow  field 
may  be  used.  Several  equations  will  be  provided  to  curve— fit 
the  pressure  distribution  over  the  various  geometrical  regions. 
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VI .  SERIES  SOLUTION  NEAR  THE  WALL 

The  introduction  of  the  transformation  variable,  \fr,  results 
in  a  system  of  equations  which  contain  a  singularity  at  the 
wall  as  a  result  of  the  vanishing  velocity,  u  -*  0.  In  Ref.  1 
the  method  was  discussed,  for  treating  the  singularity  by 
developing  series  solutions  in  powers  of  (which  satisfy 

the  differential  equacions)  for  the  variables  u,  H,  and  c.; 

l 

the  coefficients  of  the  series  were  also  shown  as  explicit 
functions  of  the  flow  properties.  Utilization  of  rhe  series 
solutions  provides  for  the  analytical  evaluation  of  the  slopes 
at  the  wall;  these  slopes  are  required  in  order  to  obtain  the 
shear  and  heat  transfer  at  the  wall. 

The  finite-difference  method  adopted  for  solution  of  the 
system  of  equations  requires  accurate  representation  for  the 
gradients  of  the  dependent  variables.  A  straightforward  point- 
wise  differencing  method  in  the  ^  direction  between  the  wall 
and  first  mesh  point  does  not  provide  sufficient  accuracy  for 
the  calculations.  The  slopes  in  this  region  are  formulated 
analytically  from  the  series  solution.  The  coefficients  of  the 
series  are  determined  by  fitting  to  the  computed  values  of  the 
variables  determined  from  the  solution  along  the  mesh  lines  of 
constant  ''J.  Terms  to  the  first  power  of  ^  will  be  retained 
co  that  only  two  coefficients  of  the  series  need  to  be  evaluated 
after  each  (incremental)  streamwise  step. 
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